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ABSTRACT 


The problem considered is the application of a De en on 
maximum principle to second order, linear differential equations of 
the form 

u!" + g(x)u! + h(x)u = f(x) for a ¢x¢ b 
with associated general boundary conditions to obtain functions z, (x) 
and zo (x) such that 
z(x) & u(x) < 2, (x) 
on [a,b]. The functions f,g and h are assumed to be bounded. We wish 
to determine the behavior of the solution u(x) on [a,b] and also to 
obtain reliable numerical estimates of u. 
The basic concepts in the theoretical background are expanded 


versions of a presentation in Protter and Weinberger [Ref. 4). 
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I. ONE-DIMENSIONAL MAXIMUM PRINCIPLES 


If u is a real-valued continuous mapping of a compact set X, then 
u(X) is compact and by the Heine-Borel Theorem, u(X) is closed and bounded. 
Thus, u is bounded. Then there exist points p and q in X such that u(q) < 
u(x) < u(p) for all x in X; that is, u attains its maximum (at p) and 
its minimum (at q). This is a specific example of a general maximum 
principle. We begin our discussion with a study of some specific maxi- 
mum principles; namely, where and under what conditions functions of a 
single real variable which satisfy differential inequalities can attain 
their maximum value. The approach parallels a similar investigation in 
Protter and Weinberger [Ref. 4]. 

Definition: For u = (x) in c? [a,b], define the linear operator L 
on (a,b) as follows: 

Ltj= u + gU 
where g(x) is any bounded function. 

Lemma 1: If u = u(x) in c? [a,b] satisfies the differential in- 
equality L [u] > O on (a,b), then u cannot attain a maximum at a point 
interior to (a,b). 

Proof: Assume, to the contrary, that u assumes a relative maximum 
atc, agc¢€b. Then u'(c) = 0 and u"(c) $ 0, hence 

Laj = uto+giouie) = uMe)<Lo 


=C 
which is a contradiction of L [u]>0 on (a,b). Hence u can attain a 





maximum only at the end points a or b. We note that the boundedness of 


g was necessary in the proof of Lemma 1 since, if lim g(x) = ©, then 
LC 
L [ul] would be undefined. The boundedness condition for g may be 
XFC 





weakened. The fact that g(x) is bounded on every closed subinterval of 


(a,b) will be sufficient for our purposes. It is important to note that 
g bounded on every closed subinterval of (a,b) does not imply that g is 
bounded on [a,b]. Consider g(x) = 1/(1-1) on (-1,1), and g(x) unbounded 
as xu 

Lemma 2: If v = v(x) in c? [a,b] satisfies the differential in- 
equality L [v] < 0 on (a,b), then v cannot attain a minimum at a point 
interior to (a,b). 

Proof: Note that -v satisfies the hypotheses of Lemma 1 since 
differential operators are linear; that is L [-v] = -L [v] > 0. 

Hence -v cannot attain a maximum interior to (a,b) which implies v cannot 
have a minimum interior to (a,b). We now have as tools two simple one- 
dimensional minimum and maximum principles. For convenience we will 
generally discuss explicitly only the maximum principles; but with 
reasoning similar to the proof of Lemma 2, we will be implicitly dis- 
cussing the related minimum principle. 

To continue our development we relax the condition of strict in 
equality on the differential inequality and now consider functions u 
such that L [u] 0. The solution u identically a constant does not 
satisfy the strict inequality in the hypotheses of Lemmas 1 and 2 since 
L [u]ž0. We now allow as a special case the constant function, which 
attains its maximum at every point in [a,b]. We will show in a Corol- 
lary to Theorem 1 that the functions u = constant are the only func- 
tions which satisfy L [u] 20 and attain their maximum interior to (a, 
b). 

In investigating the behavior of a function u which attains its 
maximum at x = a, we might imagine that u could have a graph of the form 


shown in Fig. la (and similarly if the maximum is attained at x = b). 


That is, u attains its maximum at x = a with u'(a) = 0 (or at x =b 


with u'(b) = 0). 
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(a) (b) (c) 
kipa 

We show now that this behavior is impossible if u satisfies L [u] 2 0, 
If M is the maximum of u on [a,b], then u(a) = M implies u'(a) < 0 
(Fig. lb): and u(b) = M implies u'(b) > O (Fig. lc). 

Theorem 1: Let u = u(x) in c? [a,b] be a nonconstant function 
which satisfies the differential inequality 

a! 4 gyu > O 

on (a,b) and has one-sided derivatives at a and b, and suppose further 
that g is bounded on every closed subinterval of (a,b). Then: 

(i) the maximum of u occurring at x = a and g bounded below at 


x = a imply u'(a) < 0; 


(ii) the maximum of u occurring at x = b and g bounded above at x = b 
imply u'(b) > 0. 
Proof: Suppose that u(a) = M and u(x) € M fora<xd¢b. Since 
u is nonconstant by hypothesis, there exists a point d in (a,b) such 
that u(d) < M. Define the function 
Alxa) 
zx)=e - 4 


with A a positive constant to be determined. We note that z(x) is 


positive on (a,b) and that 


~d 


7) 4 2 NY A-4) 
LJ = 2 +900200 = A € + gvAe 


Až- a) 


= A [a+ gw] e 


We choose A large enough so that L [z] > 0 on (a,d); that is, we 
select A such that it satisfies the inequality 

A >- 90 
on f[a,d] which is possible by the boundedness assumptions on g(x). 
Now form the function 


wir) = UM) + E 200) 


where E is a positive constant chosen so that it satisfies 
E < M - ud) 
2(d) 


This is possible since the numerator is positive and since u(d) ( M. 
Now j 
aia L 
on (a,d) and hence the maximum of w on [a,d] must occur at one of the 
end points by Lemma 1. We have 
w(a) = uta) + E 2a) =M 

and 

WA) < nA) + E 2td) < UA +M- ald) =M 
upon substituting the inequality for E; hence w(a) > w(d) implies 
the maximum occurs at a. Obviously, the one-sided derivative of w at 


x = a cannot be positive for then w(a) would not be a maximum. Thus 


wa) = uta)+ E zta) go 


and since z'(a) = À > Osaditefollisvs that wl(a) < O, the desired re- 
sült in (ï). 

Suppose that u(b) = M and u(x) € M fora < x 4 b. Since u is 
nonconstant by hypothesis there exists a point d in (a,b) such that 
u(d) < M. Define 


“A(L=b) 
zik) = é < 


1 


with A a positive constant to be determined. Then z(x) is positive on 
(a,b) and 


„AlE b) -~ Alx- b) 
L[e] = 2’+ m2’ = Ave ~ Age 


= A E = gin | e sosa 


We choose À large enough so that L [2] > O on (d,b); that is, we 
select A such that it satisfies the inequality | 
A > 904 
on [d,b] which is possible by the boundedness assumptions on g(x). 
Now we proceed as before and form the function w(x) = u(x) + E 2(x), and 
we reach the conclusion that w attains a maximum at x = b. The 
one-sided derivative of w at x = b cannot be negative for then w(b) 
would not be a maximum. Then 
w'(b)= utb) + E 2⁄b) 20 
and since z‘(b) = -A < 0, it follows that u'(b) > 0, the desired re- 
SUlt in (11). 
Corollary 1A: (One-Dimensional Maximum Principle) 


E 2 
Suppose u = u(x) in C [a,b] satisfies the differential inequality 


u” + g(x) u' > 0 on (a,b), with g(x) bounded on every closed subinterval 
of (a,b). If u attains the relative maximum M at an interior point c 
of (a,b), then u(x) = M. 

Proof: Assume u has a maximum at c in (a,b), then u'(c) = 0. Ap- 
plying Theorem 1 to the interval (a,c) and (c,b), we conclude that u is 
constant. This concludes the proof. We may prove Corollary lA in- 
dependently of Theorem 1 as follows: 

Assume that u(c) = M and that there exists a point d in (a,b) such 
that u(d) < M. Take d > c for convenience. Then define the function 

Å A 
Z (EAR-a) — e-a) 
where A is a positive constant to be determined. We note that: (i) 
z(x) < 0 fora < x < c, (ii) z(c) = 0, and (iii) z(x) > 0 for c< x 


< b. We have 


Liz) 


Z’ + gix) z 


A-2 TI 


A 
A(A-1)(x-a) + g(x) Alk-a) 


ij 


A-2 
= A[(A-1) +300 0-20] (4-4) 


We choose A large enough so that L [z] > O for ag x < d; i.e., we 
take A to satisfy 

(A-1) + gix) (x-a) > 0 
or equivalently 


A> -909(x-4) +14 


which is possible since g is bounded. Define the function 


wlx) ult) + E 2()) 


where E is a constant satisfying 


10 


m- uid) 


SS = 2(d) 


The existence of E is guaranteed by the hypothesis on the point d and 
the function z. Hence on (a,c) 
wit) = n) + E 24x) < M 
by (i). Also 
wld) = uld) + E 2d} < utd) +M- ua) =M. 
At the point c 
We) = ule) + E 2) = M. 
Hence w attains a maximum greater than or equal to M at an interior 
point of (a,d). But 
LW) = Ll pe ala 200 
by construction of the function z, and we contradict Lemma l. If 
d < c, then use the function 
4 A 
lix) = (a- a) — (x-a) 
which satisfies (i) z(x) > 0 on (a,c) (ii) z(c) = 0, and (iii) z(x) 
Ć 0 on (c,b). Note that the inequalities (i) and (iii) are reversed 
for this case. This concludes the proof. It would also be possible to 


use the function 


AlX—C) 
ik) = e = 
as the auxiliary function in the proof of Theorem 1, if d < c, and 
a AML O ) 
ZX) = A = 1 


if d> c. 


11 


Corollary 1B: A nonconstant u satisfying the differential in- 
equality L [u] 2 0 on (a,b) cannot have a relative maximum at an 
interior point. 

Proof: Assume, to the contrary, that u is nonconstant and has a 
relative maximum at an interior point c. Then apply Corollary IA to a 
subinterval I containing c, where I is small enough that u(c) is an 
absolute maximum on I. This yields the immediate contradiction that 
u = M. 

Corollary 1C: A nonconstant function u satisfying the differential 
inequality L [u] è O on (a,b) can have at most one relative minimum in 
the open interval (a,b). 

Proof: We note that between any minima there must be a relative 
maximum which contradicts Corollary 1B. 

The boundedness properties of g are required for Theorem 1 and 
Corollary 1A as shown in the following example. Consider the functions 
u(x) = cos x and g(x) = -cot x, which satisfy the differential equation 

U” + a(x)u'= -Cosx — “CK (- smc) 

= —lox + Cosa = O 

and here u" + g(x) u' œ 0. Conclusion (i) of Theorem 1 does not hold on 
[0, “M /2] since u'(0) = -sin(0) = 0 and u is a nonconstant function. 
On [- 1/2, 1/2] the conclusion of Corollary 1A fails to hold since 
u(0) = cos(0) = 1 is a maximum at an interior point and u is a noncon- 
stant function. We note that g(x) is unbounded in the neighborhood of 
the origin. 

Definition: A function u(x) has a horizontal point of inflection 
at x = c if u'(c) = 0 while u is strictly increasing or strictly de- 


creasing in some interval containing c. 
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Corollary 1D: A function u satisfying the differential inequality 
L [u] > O on (a,b) can have no horizontal point on inflection in (a,b). 

Proof: Assume, to the contrary, that at the point c in (a,b), u'(c) = 
0. Then on [a,c] or [c,b], u must attain its maximum on the interval at 
x = c, a contradiction to Theorem 1. 

Definition: For u = u(x) in c? [a,b] and h(x) and bounded function, 
define the linear operator (L+h) on (a,b) as follows 

(L+h)[u]= n”+ gin’ AKN 
where g and h are any bounded functions. 

In this more general setting, it is necessary to modify some of the 
preceeding Muss Lon as shown by the following examples. Consider Elie 
differential equation 

ur au = O Or e, 0] 
which has the solution u(x) = sin x. Now u attains a maximum of 1 at 
x = o. an interior point of [0,1]. 
We also Bee that the condition h(x) 4&0 is not sufficient to yield. 


r 


an unrestricted maximum principle. Consider 
an= o on [-l, // 

which has the solution u(x) = sat as and u attains a negative maximum 
of -2 at x = 0. We will prove that a nonconstant function u SSB Sn 
(L+h) [u] > O vith hG) É 0 on (a,b) cannot attain a non-negative maxi- 
mum at an interior point. 

Lemma 3: If u satisfies (Lth) [u] 2 0, with h(x) £ 0 on (a,b), 
then a nonconstant function u cannot attain a non-negative maximum at 
an interior point. 

Proof: Assume, to the contrary, that u(c) = M fora ¢ c < b and 


u(x) £ M on (a,b). Then h(c)u(c) 0, u'(c) = 0, u"(c) & 0; hence 
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te 


(L +h ra) nte) + genre) + ho) ne) 


r= 


= ue) + he) ure) é o 


a contradiction of (L+h) [u] > 0 on (a,b). 

Theorem 2: If u = u(x) in c? [a,b] satisfies the differential 
inequality 

(Lth) [n] © W+ ganu + ArU >o 

on an interval (a,b) with h(x) & 0, if g and h are bounded on every 
E sine aT, and if u attains a non-negative maximum value M at 
an interior point c, then E = M, 

Proof: Suppose that u(c) = M and u(x) € M fora<x < b and M>0. 
Then assume that there exists a point d in (a,b) such that u(d) «< M and 
for convenience take d > c. Define the function 

A (x~ a) 

| 2/x)= e = 

EEN is a positive constant to be determined. We note that: (i) 


z(x) < 0 on (a,c): (ii) z(c) = 0; and (iii) z(x) > O on (c,b). Now 


(Lrhllaj= 2/4 90)2'4 bx) 2 


2 Alt-e) Alx-e) Alx-c) 
= A e 


+ 4goje + Aitj[e — f 


= e [oče Ag) + hia] — "WD 
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We choose A large enough so that (Lth) [z] » O on (a,d), that is, 


so that A satisfies the inequality 


X _A(x-€) 
[a + A9) +hto]|- huje > 0, 
Since SPA and h & 0 on (a,b), we need only to choose A such that 


A*+ Agqit) +hix) DO 


which is possible since g and h are bounded. Now form the function 


Wit) = Uk) + E 2K) 
where E is a positive constant so that 


M- nla 
ES ma) 
Hence on (a,c) 


Wht) = HALE DL) < M 
since z(x) < 0 on (a,c). And by the definition of E, 


wid) = Uld) + M— uid) = M 


Now 
W) = HCO) + E Zle) = M 
since z(c) = 0; thus, we may conclude that w has a maximum greater than 
or equal to M which is attained interior to (a,d). But 
(Ltn) lw] = (Leh)lap + E (Lehllz) > 0 
on (a,d) and we have a contradiction to Lemma 3 since w(x) is noncon- 
stant. 
If d € c, construct the function 
2/1) s ika 1 
which has the properties: (i) z(x) > 0 o (a,c), (ii) z(c) = O, (iii) 


z (x) < O on (c,b). Now 


15 


(L +h) [2] = | 2+ 91k) 2° + Alki 2 


~ 


2- Alx-¢) ~ AX?) ~Alx-e ) 
pe“ gaga” ahaj <a) 


e ee [A Aq (x) t hraj — Atk) 


We choose A large enough so that (Lth) [z] > 0, that is so that the 
following inequality is satisfied: 
Z Alk-ca) 

[ K- Agw + huf- Aix) e > 0 
~ Alx-c) 


since @ > 0 on (a,b). And since h & 0, choose A such that 
A Agik)+hix) DO 
which is possible since g and h are bounded. We note that the condition 
that A satisfies 
A+- Alg + hix) > o 
will suffice for both d € c and d > c. The remainder of the proof is 
exactly the same as before with c < d. 
Corollary 2A: If u satisfies (Lth) [u] 2 0 on (a,b) with 
h & O and h(x) SÉ O, and u assumes a non-negative maximum value at an 


interior point, then u(x) = 0. 
Proof: Let M & 0 be the maximum attained by u on [a,b]. By 


Theorem 2, u(x) Æ M and hence 


(+h) [M] = hj M 2 O 
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But h(x) É O and h(x) £ O implies M = 0. 

Theorem 3: Suppose that u = u(x) in c? [a,b] satisfies the dif- 
ferential inequality 

(Lehjln] = Ut GAIN + hwn >o 

on (a,b) and has one-sided derivatives at a and b, that h & O and that 
g and h are bounded on every closed subinterval of (a,b). Then: 

(i) the non-negative maximum of u occurring at x = a and the func- 
tion g(x) + (x-a)h(x) bounded from below at x = a imply that u'(a) < 0; 

(ii) the non-negative maximum of u occurring at x = b and the func- 
tion g(x) - (b-x)h(x) bounded from above at x = b imply that u'(b) > 0. 

Proof: We may extend the proof of Theorem 1 by noting that if the 


maximum of uc occurs at x = a, then 


AK-A 
Aix- A) 
(L+th)[e - {|= e 


) = 9 1-4) 
[4% AG ix) + ha)(1-e y 


Alx- 


} 
2 e Jače Ag) + A (x-a) hra) 


H / 
since@ en for t > O and h(x) £ 0 by assumption. If g(x) + 


(x-a)h(x) is bounded from below at x =a, then we may choose A such 


' A» -[300)+ 1-4 ht] 


If the maximum of u occurs at x = b, then as above 


Alk- —Alx-b) = 
(cone? Je e” [pr-dgyt baste” Y] 


— AlX-&) 
Deo TA Agu) + A(b-x) hia] 


Ly 


Hence if g(x) - (b-x)h(x) is bounded above at x = b, then we may choose 
A such that 

A> 9) ~ (b-x) hik) 
The remainder of the proof is analogous to Theorem 1. 

Corollary 3A: If u z O satisfies (L+h) [u] 2 0 on (a,b) with 
h 0, if u is in car and if u(a) & 0, u(b) & 0, then u(x) < 0 
on (a,b). 

Proof: Assume that there exists a point c on (a,b) such that u(c) 

> 0. Then u(a) £ O and u(b) & O imply that u attains a non-nega- 
tive maximum at some point d interior to (a,b),say u(d) = M. Then, by 
Theorem 2, u(x) Æ M BO, a contradiction. 

A useful technique in determining the existence of a maximum (or 
minimum) of a solution to a differential equation at an interior point, 
independent of the preceeding theory, is a simple direct argument. 

Consider the problem 

nr e” = —x on (ot) 
which does not have the format of the problems we have so far considered 
due to the exponential term an 

We note that 

fu) = w+ ze = -X <0 
on (0,1). Hence, u cannot attain a minimum in (0,1); otherwise, if 
there exists a point c in (0,1) such that u(c) = M and u(x) 2 M for 


x s c; then u''(c) > 0 and we have 
Zle) 
fm) = uvlejtre >o 


on (0,1), a contradiction. 


As another example of the preceeding, consider 


A"—~2 G09 (u’) = 1 
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whose solution u cannot attain a local maximum in the interior of any 
interval. Let I be an interval with subinterval [a,b]. If there exists 


c in (a,b) such that u attains a maximum in [a,b] at c, then 


u (0)e0 and “ve ) = O 
Hence 
O> U") = 14 2c05/uto] = 14+ 2005 (0) = 3, 
an impossibility. 
We may even apply the previous results to boundary value problems. 
Consider B xy 
H+ Cee ee i 
for 0 x 1, with u(0) = u(l) = O. Now 
u'selu'<o 
on (0,1); hence u cannot attain a minimum interior to (0,1) as in the 
first example. The boundary conditions u(0) = u(l) = O imply that u'(0) 
> 0 and u'(1) £ 0; otherwise u would attain a minimum interior to 
(0,1). If u'(0) = O and u(0) = 0, then from the differential equation, 
u"(0) = -1, The boundary condition u(l) = O forces u to attain a mini- 
mum interior to (0,1), a contradiction. We obtain a similar result if 
we assume u(l) = u'(1) = 0, since u(0) = O. Thus u' (0) > O and u' (1) 


Lo. 
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II. A GENERALIZED ONE-DIMENSIONAL MAXIMUM PRINCIPLE 


If we drop the restriction that the function h must be non-positive 
we may then obtain a maximum principle for (L+h) [u] > O by introduc- 
ing a function w, with certain desirable properties, such that the 
quotient u/w satisfies the previous maximum principles. 

Theorem 4: Suppose u = u(x) in c? [a,b] satisfies the differential 
inequality (L+h) [u] 2 0 on (a,b) with h(x) bounded and g(x) bounded 
below; then for any sufficiently small subinterval [a',b'], there exists 
a function w in c? [a,b] such that: 

(i) w> Oon [a',b'] 

(ii) (Lth) [w] £ 0O on (a',b'). 

Then the function u/w satisfies the maximum principles of Theorems 2 and 
3 on (a',b'). 
Proof: Consider the function 
2. 
wix)a 4- Blx- 4º) 
fora «a' < x < b' < b, where B is a positive constant to be 
determined. When B has been determined, we will suppose that (b'-a!) 


is such that 
2 
(1) B(b'-a') < L 


so that w > O on [a',b']. Now 


(LZ +h) [w] = Ww" 91x) w’ + hix) w 


= (-28) + qumf-2e aj + Alt) J. e Bla) 
= -28 / 1+ dra) Gia) + i lia’) hn) + Aik) 
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Since g and h are bounded below, choose constants G and H such that 


g > Gand h 2 H. We now make a restriction on (b',a') so that 
E i 3 
O 1+hx-4)Q + 2 (4-4) #>0 


fora! £ x € b! and choose B so that (Lth) [w] & O on (a! ,b'), 


dese; hik) 
Dee pi ii 
£ 1+ a) GA le-a y 
a permissible choice since h is bounded above. Hence if the subinter- 
val [a',b'] is sufficiently small, the function w with properties (i) 
and (ii) may always be constructed. 


Define the function v = u/w on [a',b']; then 


(hath) [up = (lerh Ere] = (ru) e gen tru) + hli vw 


= /ru"t2vuruv"]e gut) (rare) shi) ru 


Now dividing by the positive function w, we obtain 


/ yy! Pd A > 
ge + [2% + guy iP 5 (hth) lw] ¥ = 0. 
Thus, v satisfies the differential inequality 


T+ Imr + Alu 20 
on (a',b') where 


g, ix) = 2M + HK) 5 hiti~ Zh) lw) | 


a 


By the construction of w, w ín še [a',b']; thence g and h are bounded 
on [a',b']. Also h is non-positive by conditions (i) and (ii). Thus 
v satisfies the hypotheses for Theorems 2 and 3 on subinterval [a',b'], 


which concludes the proof. 
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We now see that a function u which satisfies (L+h) [u] 2 0 has 
certain restrictions on its zeros on the interval under consideration. 

Corollary 4A: If the function u satisfies Theorem 4 on (a',b') 
and u(x) # O, then u can have at most two zeros (between which u is 
negative) on (a',b'). 

Proof: We may take the interval [a',b'] such that u(a') = u(b') = 


O and u(x) + O fora! < x < b. Assume that u(x) > O for a! & xd b. 


Then v(x) = u(x)/w(x) > O fora! < x ¢ b' and v(a!) = v(b!') = 0. 
Hence v has a non-negative maximum in (a',b') which is a contradiction 

to Theorem 4. Thus, u(x) < 0 on (a',b'). Also, there could not be more 
than two zeros; otherwise, if there exists a point c in (a',b') such that 
u(c) = 0, then u(x) < _0 on (a',b') forces u'(c) = 0. But then u must 

be the trivial solution, a contradiction. 

Note that the zeros of u must be isolated; i.e., if u(a) = 0, then 
u cannot vanish on some finite interval to the right of a. 

Corollary 4B: If u is a solution to the differential equation (Lth) 
[u] = O, then u can have at most one zero in any interval (a',b') where 
Theorem 4 holds. 

Proof: Note that -u also satisfies (Lth) [-u] = 0. By Corollary 
LA, u and -u can each have at most two zeros on (a',b'), between which 
both u and -u must both be negative, an impossibility. Thus, u has at 
most one zero on (a!,b!). 

We may conclude that the M undedneik hypothesis in Theorem 4 on the 
function g is essential by noting the following example. The function 
u(x) = x satisfies u" +g(x)u! = 0 vith g(x) = -2/x, and g is unbounded 
in the neighborhood of zero. Consider the interval of the form (a,b) 


with b = O. Then u attains a non-negative maximum of zero at x = b and 
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u'(b) = 0, contradicting conclusion (ii) of Theorem 3, hence contra- 
dicting Theorem 4. 
Definition: Let r(x) be a solution of the differential equation 
t+ GM + ht =O 
with g and h bounded functions. If r z O and r(a) = 0, the r must be 
non-zero on some finite interval to the right of a by Corollary 4B. If 
r has another zero on (a, e ), denote the first such zero by a* and call 
it the conjugate point of a. (If r £ O on (a, “O ), then set at = 25), 
Thus r has one sign on (a,a*) and for convenience we may take r(x) 
2 0 fora < x< a*, since -r satisfies the same differential equation. 
We now establish a necessary and sufficient condition for the existence 
of the function w of Theorem 4 in terms of a*, 

Corolllary 4C: If a* is the conjugate point of a, there exists a 
function w > 0 such that Theorem 4 holds on the interval [a,b] if and 
only inb € at. 

Proof: Assume that b > a*, Ifw > 0 on [a,a*], then the func- 
tion u/w vanishes at a and a* since r(a) = r(a*) = 0, and u/w is posi- 
tive on (a,a*). Hence u/w attains a non-negative maximum in (a,a*) and 
u/w is non-constant. then w cannot satisfy (Lth) [w] < 0 on (a,a*) with- 
out contradicting Theorem 4. Hence by a contrapositive argument, we 
have the forward implication that the existence of an acceptable w(x) 
implies b «< a*. 

To get the reverse implication, let b in (a,a*). Note that r(x) is 


bounded from below by a positive number m = m(c) on any subinterval [c,b] 


witha < c. Let 4tx—a) 
Wh tac - Ee 


where A is a positive constant to be determined. Consider the function 
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wiz) Y) + E 21%) 


Now choose E small enough so that w > 0 on [a,b]. Now 


(1th) [w) = (bb) + sluze") 


> ALA) A (x-a) Alx-4) 
= E fife ) +g (-Ae  )+ hl (2-e iy 


=— Ee SAS A gua + ha) ji - že“ a 


We want to choose A such that (Lth) [w] & 0 on [a,b] which can be done 
as follows: 
Alt-a) A (ia) 


Now -Ee < 0 since E and e are positive, We make 


two restrictions on A: 





- Ali-a) > /n 2 
(i) 1 - 2e PROMBIVE D4 A 
a 
Since h is bounded, we may also choose A large enough so that 
— A (X-A) 


GD A+ Agik) + hw |4- Ze po. 


Then we have (Lth) [w] & 0, and w is a function for which Theorem 4 holds. 
If h(x) is unbounded or g(x) is not bounded below there may be no 
interval [a',b'] for which Theorem 4 holds. Consider the function 
nik) = E sm (É) 
which satisfies 
nsx H =O 
Now as x—> 0+, h(x) = x—» + Se“, There can be no function w > 0 
with the property that u/w satisfies the maximum principles in any of 
the intervals [0, 1/n9” ] (n = 1,2,...) since u(l/n4”) = 0 and 1/nùy>0 
as n—> So , Thus, for all ¢ > 0, the interval (0, & ) contains a 


zero of u and Corollary 4A denies the existence of an appropriate w(x). 
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As in the examples in the previous section, we use a direct argument 


to conclude that no solution of 
yf as 
A re 
can attain a minimum value greater than 1 or a maximum value less than 1. 


Assume there exists a point c such that: 


(i) u(c) = maximum; hence u"(c) < 0 and we have 
Hie) 
e — e =~ y") 26 


Ue) 
which implies e  >e, and we have the maximum value of u is u(c) > Infe) 


=]; 
(ii) u(c) = minimum; hence u"(c) 2 O and we have 
u (e) A 
e -a=u(c) £0 
ute) 
which implies that e <= e and we have the minimum value of u is 


u(c) £ ln(e) = 1. 


III. INITIAL VALUE PROBLEM 


Consider the solution to the differential equation 
4 eet h one 
œ u'rgwu+ hwu = +0 


which satisfies the initial conditions 
; I 
oO rUMD=% 5 Ua= G 


where the functions f, g and h are defined on the interval (a,b) with 


g and h bounded and and LA given constants, 


Y 
1 

The existence of a solution follows from the classical theory of 
ordinary differential equations [Ref. 3]. We may prove the uniqueness 
of such a solution independent of the classical theory using Theorem 4. 

Theorem 5: Suppose the u, (x) and u, (x) are solutions of (1) which 
satisfy the same initial conditions (2) in an interval (a,b). Then 
u = u, in (a,b). 

Proof: Define u(x) = u, (x) “uy (x). Then u satisfies 

U"+ 9X) U'+ HA)U =O 
with initial conditions 
Alah= Ua) =O 

Now assume that u + O in (a,b). Then by Theorem 4 there exists a 
positive function w on an interval (a,a+€ ) such that the maximum of 
u/w in (a,at& ) must occur at one of the endpoints. We note that 
is a positive constant depending only on the bounds of g and h. I.e., 
(a+ €) - € = € must satisfy conditions (l) and (2) in the proof of 


Theorem 4. But -u also satisfies the same equation with the same initial 


conditions, hence, also by Theorem 4, -u/w must also attain its maximum 
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at one of the endpoints, a or at. Hence u/w must attain either a maxi- 


mum or a minimum at a. But at x =a 


/ 


(=) u'w- Uw! 


u = a = 


A Ww + Lea. 
since u(a) = u'(a) = 0. Now the function u/w satisfies Theorem 3; hence 
u/w must be identically a constant, and from the initial condition, u =0 
on [a,a+€ ]. In particular 

| 7((01€) = u'(at£) =0 
Hence, by the analogous argument, u= 0 on [até , at2& J with £ un- 
changed since the choice of € is independent of u. Now (b-a)< = ; 
hence, after a finite number of steps, we may ona that u = 0 on 
[a,b]. Thus UE u, 


Consider the problem 


u = O 


/ 
U 
K x 


with the initial conditions 


yio) = U0) = O 
We note that uj = 0 satisfies the equations trivially and for u, = e 
7, — + Ne + (2x) =e 
But u, z u, on (0,1). Since g(x)—> =© as x—> 0, we may conclude 
that the boundedness of g is necessary in Theorem 5. 
We can easily find a function h(x) such that the equation 
Doe RA 
has two solutions satisfying 


4lo) = 41o) 


1 


O 


Obviously u, =Œ 0 is a solution. We must then pick h(x) unbounded in 
the neighborhood of the origin; otherwise Theorem 5 would imply that the 
trivial solution was the only solution. Choose h(x) = “6 fe" then u(x) 


= x is an appropriate solution satisfying the initial conditions. 


Hence we see that the boundedness hypothesis on h is also necessary in 


Theorem 5. 


An interesting application of the Mean Value Theorem arises in 
the initial value problem 


(3) H+ a gy for x>o 


satisfying 
(4) wos=l and y7oj= ©. 


Assume that u) and u, are solutions of the differential equation (3) 


satisfying the same initial conditions (4). Then the function 
HIK) =E U,K) DUE) 


satisfies the homogeneous differential equation 


(H HA) + Fe. e 2) = Ø, 


By use of the Mean Value Theorem 
d, uz Hz + Olid) 


e-e = (tz) € /o0<oe| 


and the differential equation becomes 


e+ OCU; -Up) 


(4,'- U,") +,-h, ) e = o 
with 
dio) = 14/60) =0 
h ÒH- A) 
Since h(x) = -e <= 0 and is bounded on any finite sub- 


interval of (0, oo 5, we may apply Theorem 5. The trivial solution 
u Æ O satisfies the above initial value problem, thus u is unique 


i = 2 Hence we have the solution of (3) satisfying (4). 
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IV. BOUNDARY VALUE PROBLEM 
We will first be concerned with boundary value problems of the 
type 


Do onegmaehnu = Fix) 


on the interval (a,b) subject to the boundary conditions 
gy wal= 4 5 UW Kn. 


We may derive a uniqueness theorem by the use of Theorem 4 for solutions 
of (1) which satisfy (2). The sign of the function h(x) has great signi- 


ficance as seen in the simple example: 


ha Tr 


which has the solutions uy = 0 and u, = sin x on the interval [0,1]. 
Both Wy and u, satisfy the boundary conditions u(0) =u( MK) = 0. 


Theorem 6: Suppose that u, (x) and u, (x) are solutions of (1) 
which satisfy the boundary conditions (2), If h(x) € 0 on (a,b), then 


u = 


1 2: 
Proof: Set u(x) = u, (x) - u, (x); then u satisfies 
ata gu) ut hu = 
with the boundary conditions 


Hla)= ulêl = O 


By Corollary 2A, u(x) £ 0 on (a,b). By analogous reasoning -u satis- 
fies the above equations; hence «u = 0 on (a,b). Thus u = O on 
(a,b). 


In order to incorporate a larger class of boundary value problems 


we consider more general boundary conditions which have (2) as a speciai 
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case, We consider solutions of (1) which satisfy the boundary condi- 


tions 


Ra) loso + ula) dun © = TF, 


(3) | 
(bl) o0sk + ud me = |F 


where O, b E Sa 1 and (AR are prescribed constants with 0 £ 0 < ‘ 


x 
ost ; We have (2) when@= Æ =". 


To derive the general boundary conditions (3) consider the boundary 


conditions 


A, (~n'ta)) Pm, ula) = g) 
Aa U) + B UO 


where A B.» C,» (i = 1,2) are given constants. We may normalize the 


LE 
above, i.e., B é 
A / 
VA) PRE MAR) po a 
V A pe 6,* ( / VA? 6, * 7 Y Afr 6,2 7, 


Az a 


Se E W Ch é 
y ARE ME, e ) a Vag , wl k) = VE apa < a 


Then we define the angles © and : so that 


and 


Cos FAT, SU | E. 
= J = 
, & At 8,2 2 no A +6)? 
an 
A, B 


4 


Cosp = VAZA E , Sin & 
4 
Then 04 6474 and 0 £ BE . 


Vaz + #2 


Theorem 7: Suppose that u, (x) and u, (x) are solutions of (1) 
which satisfy boundary conditions (3). If h(x) = 0 on (a,b), then 
uy = u, unless h Es O q A £ = 0; in which case u, and u, may differ 
by a constant. 


Proof: As before set u = uj 7 Uys SO that u satisfies 


W"+9Q)0' Ali u = O 
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and 


4) -UVA) CO + HM) IME 


it 
Sy 


6) ale) est + Hl) Sind 


if 


O 


Consider u = M, a non-zero constant. Then (4) and (5) imply 
Msme= MsngG=0 

Now M 0; hence sin& = sing = 0 and o=¢4 = O. By Corollary 2A, 

M #0 if and only if h= O. Thus, if uj and u, differ by a constant, 

then h = 0 and & = $ = 0. 

Conversely, assume © =f = Q andh = O; then u = M satisfies (4) 
and (5) trivially and 

U"+ Glt)U' + Atk) Hn EO, 
Hence, uy and u, differ by a constant. 

Next assume u is a non-constant function which is positive at some 
point in [a,b]. By Theorem 2, u attains its positive maximum at a or 
b. Suppose the maximum occurs at x = a. Then by Theorem 3, u'(a) < 0. 
Now u(a) > O and ofo4et g oe fs x * hence 

~H'CA) C050 +7 /a) wne >O, 

since cos & and sin @ cannot vanish simultaneously and are positive 
on the interval under consideration. This inequality is a contra- 
diction of (4). If the maximum occurs at x = b, we have u'(b) > 0, 
u(b) > 0, and 

7016) tos f + Wb) 3nd So 
a contradiction of (5). Hence u is a non-constant solution which can 
never be positive. Analogously -u is also a solution which satisfies 


(4) and (5); and we may conclude that u can never be negative. Thus u 


SF 0 on [a,b}. 
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Now let the only assumption on the function h(x) be that it is bound- 
ed on the interval under consideration; i.e., we remove the restriction 
that h(x) be non-positive. We will now prove a uniqueness theorem, but 
as the following example shows, the interval considered is of vital im- 
portance. The simple equation 

n"+H =O 
with u(a) = u(b) = 0 has only the trivial solution for b-a < TY ; and 
for b-a = 9” uniqueness failed as was shown at the beginning of this 
section. 

Theorem 8: Suppose that u, (x) and u, (x) are solutions of (1) 
which satisfy the same boundary conditions (2) with h(x) & 0 on [a,b]. 
If b < a*, where a* is the conjugate point of a, then u = u. 


Proof: Let u = uj = u95 then u satisfies the equation 


n'+ Hu + hinu =0 

with 
Wa) = Hb) =O. 

Since b < a* by hypothesis, Corollary 4C implies that there exists a 
function w > 0 on [a,b]. Hence by Theorem 4 the function v = u/w satis- 
fies the maximum principles of Theorems 2 and 3 and their Corollaries, 
Now v(a) = v(b) = 0 since u(a) = u(b) = 0, hence v(x) & O on [a,b]; 
otherwise, v would attain a positive maximum in (a,b), a contradiction 
E Theorem 2 and Corollary 2A. Since w is a positive function, u(x) & 0 
on [a,b]. We may apply the same argument to -u since it satisfies the 
same differential equation and boundary conditions, and we conclude that 


-u É 0. Thus u =O. 
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V. APPROXIMATION IN BOUNDARY VALUES PROBLEMS 


We are concerned with the same boundary value problem as in Section 


IV. Thus, we seek a solution of the differential equation 
=. / 
ay (Lth) [H] S me gn hiiu = fit} 
on the interval (a,b) which satisfies the boundary conditions 


o ual= Hd ; n= ha. 


Even in many of the simplest problems an explicit expression of 
the solution is impossible to determine. We must then resort to series 
solution techniques or other essentially numerical methods, Obviously 
an approximation solution technique which yields an explicit bound for 
the error compared to the exact solution is desirable. It is possible 
to find such a bound by finding explicit approximation functions which 
provide an upper and lower bound for the actual values of the solution. 

In this section we wili develop necessary conditions for the exis- 
tence of such approximation functions, and in Section VI we will in- 
vestigate more closely the approximation functions and related results. 

We assume, as before, that the functions f, g and h are bounded on 
the interval under consideration. With these assumptions it is possible 
to determine bounds for a solution u of (1) and (2) using the maximum 
principle in Theorem 3 without any further knowledge about u, 

Theorem 9: Suppose that u is a solution of (1) satisfying boundary 


conditions (2), and that h(x) 0. If the function z, (x) satisfies: 


o (Lth) [z] £ Fix) on ah) 
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(4) 2/4) > 7, and 2(b) > E 
and the function z, (x) satisfies: 


6) (Ath)[2,7 > Fr) on (4 b) 
6) alimen Ad Alb) EK, 


then 


22164) € HM) € 218) 
fora €x& b. 


Proof: Consider the function 

WO E AMI- PA) 
Then v, satisfies 

CUPA IIR T 20 
on (a,b) with 

(0) co Aand Wh) £0. 
Apply Corollary 3A; v(x) < 0 on [a,b]. Thus, 
Hy) e DK). 

Similarly the function 


Ta (k) = 24) — (x) 


(Lt) [vz] =O 
on (a,b) with 

mta) co and n, h) £0. 
Thus Vo (x) = 0 on [a,b] and we have 


2,14) £ UK) £2,1K) 


satisfies 


on [a,b]. 


A more general two-point boundary value problem is 
oy (hth) [Ww] = +930 hau = iš) 


on (a,b) where the solution u satisfies 


- Ha) coso + UG) Ine =F, 


pi 2 (4) Cosg + uib) Sind = Za 
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The quantities © and + are pre-assigned constants wuch that 
ose stag and begeg 
Continuing to assume h(x) £ O, we may extend the proof of Theorem 9 
to get the following. 
Theorem 10: Suppose u is a solution of (7) satisfying boundary condi- 
tions (8). Suppose also that h(x) & O and @ , A are in [0,^Y /2] 
and that not all the equalities h T 0, 6 =0,, f = 0 hold. If the func- 


tion z, (x) satisfies: 


o) (Lh) [e] £ Fit) en (44) 
- 2 (4) Cosb + 2,f4) Sene 27, 


(10) 


2z bless + zgh) shg ZR 
and the function Z, (x) satisfies: 

an /4+h) iz) > Ft) 04 @b) 
-3 a) DIO + Ba) me ET 
2%) arb + lb) Sind £%, 


(12) 


then 


AA S MIs 2,/X) 
Proof: Set the function 


(44h) m] 


v. Z u -z.. Then v, satisfies 
1 1 1 
> 


O 


and 
~ 2) A) 056 #40) “m6 £0 
nC ish + Xb) Bind £0. 


TE vi is always positive, then Theorem 2 states that the positive maxi- 
/ 
mum of v} occurs at either a or b. If it occurs at a, then v; (a) = 0 


and v, (a) EO. But 


_ Va) 036 + Va) Smo £0 


can occur only if sin @ = 0, i.e., © =0 and v'(a) =0. Then Theorem 
3 states that v, (x) is identically a positive constant. Thus h(x) = 0 
by Corollary 2A. Similarly Vi cannot attain a positive maximum at b 
unless g = O and h = 0, since 

xh) cosh + bl sing £0. 
We may conclude that unless © = 0, f = 0 andh Ž 0, v(x) <= 0o 
Carb): ivê 

nk) € VE). 


By an analogous argument, the function Vo = Z =u satisfies 


[Lth) [1] >0 en (uk) 


and 


- va) ie + Gl) Imo e0 


GL) EE + YO)IWG €0. 


If either h g O or not both & and É are zero, then we may conclude 
Z (x) S$ yix) E 2X) 
on [a,b]. 
Remark: If the special case h @ 0, © = f =0 holds, then from 
the proof of Theorem 10 we have that the functions vi and v, are both 


identically positive constants, say 
az= M, >0 and ghz M, >o. 
Then 
20t) LR tM, = A= Bia M $al; 
i.e., we have 
2 x) < nik) < 244) 


and Zi and Zo differ by a positive constant. 
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We drop the restrictions that h £ 0 and that &, É are in [0, /2]. 
From the discussion on the background of the generalized Sturm-Liouville 
boundary conditions we may have just as easily taken @, fs in (-%4 /2.1/2] 
with no loss of generality. 

Theorem 1l: Suppose that u is a solution of (7) satisfying boundary 
condition (8) with - M /2 < 8 E Ty? AY /2 < pf < M /2. Let 
the functions z, and z, satisfy the inequalities (9), (10) and (11), (12) 


l 2 


respectively as in Theorem 10. Then the bounds 
%, AS Ulk) E BK) 
hold in (a,b) if and only if there exists a positive function w on [a,b] 


satisfying 


a3) (bth)lw]< o em (n bj 
-wla Cose + la) *m6 >20 
(h) asp + wib) ShÊ 320 


in such a way that not all the inequalities in (13) and (14) are equalities. 


(14) 


Proof: Suppose there exists a positive function w on [a,b] which 
satisfies the inequalities (13) and (14). Then define the function v = 


u/w. As in Theorem 4 


(hta) = (Lt) [rw] 
= ve |Z a + 34) ari [tb (3h Jud fr x = 


SPA Tee 


where G = (2w'/w) + g and H = 1/w(L+h) w. Note that H, the coefficient 
of v, is non-positive. The functions v satisfies the boundary condi- 


tions? 
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1 


ty YA) WA) YA) wa) COSG + V4) 4A) SMe 


x, 
Va 


il 


[orca wo) + vlo) wlec + Wo) wulh) Im q 
which we may rearrange to form 
~Y¥'(A) WIA) GOS © + ya j-w (4) Wot wla) m e] =F, 


(16) 
bato) oast Vo |W Couga ulesing] = 


Dividing through by w(a)cos & and w(b)cos » then multiplying by the 
non-negative quantities cos © and cos #4 in the first and second condi- 


tions of (16) respectively, we have 





/ Ce (-ult) @e@ + Wlh) Fme) 4, C39 
VA CEE + UA) | ———— /= 





Wa) COs ula) (os 8 
(17) 
_ «Me | Fy G, ap 
(6) au + vio) Cos F Cub) d sub) If) | da 
Wb) Cos ilhi cosg 


where the angles & and A are defined by 


CosE (wa) Case + UMA) Sim + ) 


din O = 
wti) CoJo 
a cosé (wh) Usd rub) ng) 


Wb) cos 
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or equivalently 
a — u'(a) Cose + wa) sine 
Zan o = —— i a 


Wla/ CASO 
wilbpooss ~ wl6/ Sm Kk 
Be) Cos É 

By the hypotheses on & af and w, both tan 


(18) 


Ton ge = 
and tan E are non- 


o 
negative. By (14) we may always take & and op such that 


tp 


es6<% and Ocf oy, 


Note: If @ = V/2, take B = Tv /2; if A= /2 take $ = Ù /2. 


Hence we have 


vale + VA me 


H 


%, 


(19) 


{i 


yb] wsh + V) sgimh = h 
with © , $ in [0, I /2] and 2, = (7 cos 6) /w(a)cos® and 4, = 
(M cos & ) /w(b)cos fp 


If the functions Zi and 2. satisfy the conditions (9), (10) and 
(11), (12) respectively; then the functions z,/w and z,/w satisfy the 


analogous conditions with respect to equation (15). Hence by Theorem 10 


the following inequality holds 

Z lh) MH ZA, zi g, Cx) 

AK) = ‘Wik ~~ Wk) 
unless © = K = O and H E O. From the proof of Theorem 10 this situa- 
tion holds when the inequalities (9), (10) and (11), (12) are equations. 


Hence, if there is a positive function w(x) satisfying (13) and (14) but 


such that not all are equations, we may conclude 


(20) BM) E UK) CLA) 
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Note: If w satisfies the boundary value problem 
(l+b)l4] =0. On (ab) 


-wral Wo + wla) GmHO =O 


with) cost e wlb) sme = 0 


then let & be any real number and u a solution of (7) which satisfies 


(8). Then 


(Leh) [as ew]= (4rhjlii] + § (hth) Pate 


-ua Eura] Ose + | usa)+ § W| Sme = F 


[ules § WO) cos + | nce p Euo) sad = XK 


Thus, if w satisfies (13) and (14) as equations, then any multiple of 
w, when added to a solution of (7) and (8), will yield another distinct 
solution. Hence, if at least one solution u exists, then there are 
many. 


Now consider the boundary value problem 
ay (lth)lu]=o0 m (ab) 


with : 
~ WE) (0S 6 t+ LE) sime = 1 
(22) 


y's) Cosh+ wh) sm db = 1 
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By applying the classical existence theorems, let w be a solution of (21) 
satisfying (22). Noting that the function w satisfies analogous condi- 
tions to the function u in Theorem 10, we have | 

2, (x) £ Wie) & 2, (k) 
for appropriately chosen bounding functions z| (x) and Z, (x). Also 
z, (x) = O satisfies condition (11), (12); hence the function w(x) is 
non-negative on [a,b]. We now show that w(x) is actually positive on 
lamb). 

If w is zero at an interior point, then w' is also zero there, 
otherwise w would have some negative values. But the identically zero 
function satisfies the same initial conditions as w(x) and by Theorem 
5 for uniqueness of initial value problems, w E 0 on [a,b] which con- 
tradicts (22). Hence w cannot vanish at an interior point. 

If w vanishes at an endpoint, say at a, then from the first condi- 
tion of (22) we have 

WA) ose =-1. 
Now w non-negative implies that w'(a) 2 0, a contradiction since cos O 
> O. Thus w(a) > O. Similarly w(b) > O, and we have w > 0 
on [a,b]. 

As with Theorem 1 and Corollary 1, it is often possible to prove a 
more general theorem and by a judicious choice of parameters obtain a re- 
lated, more specific result. For example, if h £ 0, 0 £ 6 < ‘ie /2. 

0 < pe “i /2; then the function w = 1 satisfies (13) and (14) and 
Theorem 10 follows immediately from Theorem 1l. 


The fact that the function w does not appear in the inequality 


Cox) E UK) & 2,1%) 
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prompts us to try to obtain a necessary and sufficient condition that 
guarantees that the functions zy and 2, form upper and lower bounds for 
u independent of a function w. 

Theorem 12: Suppose that the functions Zz, (x) and z, (x) satisfy in- 
equalities (9), (10) and (11), (12) respectively, in such a way that 
equality does not hold in all the conditions. Let g(x) be bounded be- 


low on each interval [a,c] and bounded above on each interval [c,b] 


witha < c < b. Let u(x) be a solution of (7), (8). Then the bounds 
(23) 2,1) é Alx) £ Z, (x) 


hold if any only if z, (x) £ z (x) for a< x < D. 
Proof: Trivially, if (23) holds, then z, (x) = z (x). Hence assume 
Zz, 7 25 is non-negative and we must show that (23) holds. Put 
IX) = (E) — (X) 
If q > 0 on [a,b], then q satisfies the conditions for the positive 


function w in Theorem 11; i.e., 
(24) (4th) [9 ] = (ith ) f2,) - (L+h)[2,] Lo on aih) 


since 2| and 2, satisfy (9) and (11) respectively. 


~[ 2,4) — 2 af ose + [2 täi- Gla] sme So 


[2/7] - 2 (b)/Cosd + / 2,0) - 2,04] sine 20 


since 2, and Z5 satisfy (10) and (12) respectively. Hence we may con- 


(25) 


clude by Theorem 11 that 

ZAK) & WK) e 2x), 
Thus, we need only consider the case where q has a zero in [a,b]. We 
will show that q cannot vanish at an interior point and that if q 
vanishes at one or both endpoints then we may find suitable positive 


functions w(x) such that Theorem 11 will hold, 


4? 


Suppose that q(c) = O for c in (a,b). Then q'(c) = O since q > 0 
on [a,b], i.e., q has a local minimum at c. We may conclude that q = 0 
on [a,b] by the uniqueness of solutions to the initial value problem. 

But then equality holds in (24) and (25), a contradiction. 

We are left with q > O interior to [a,b] and q vanishing at an end- 
point, say at x = a. Then q'(a) > O by above. Thus & = IRA 2, other- 
wise the first condition of (25) would be violated. Similarly if q(b) = 0, 
then q' (b) < 0 and ob = Y /2. 

We now consider the case where q(b) = 0 and q(a) > O, and proceed 
to find a positive function w such that Theorem 11 will hold. From the 
preceeding paragraph B = TM /2. Consider the initial value problem 

(4th) [VY] = O 

(26) , i 
Pa) = Ce and Via) =sne 
which has a solution by classical existence Theorems for differential 
equations. Let r be the desired (unique) solution. We note that for 

Q < dao. r is (strictly) positive at a, and if O = dia /2, r(a) 
= O and r is positive near a. | 


Now define the function v = r/qfora<x< b. Also note that 


a)? = = 
v(a) T] QO 


rija O VA) - Za) Wa) 


Fla) 
70) Fine — Pla) Cole 
Za) 


by the first condition of (25). Using the facts that (lth) [r] = 0 


and 





q 
wan 


and (Lth) [q] O, we may conclude that v satisfies 


(LtH)fuJ=u4av's+ A =O 
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with G = (2q'/q) + g and H = 1/q(Lth) [q] £ O, as in the proof of 
Theorem 11. The coefficients G and H satisfy the restrictions in Theorems 
2 and 3 on any subinterval of the form [a,c]. Hence either: (i) v(c) 
> v(a) and v' (c) > 0, or (ii) v(x) = v(a) for x in [a,c]. In either 

case r(x) > O for x in [a,b) since both v and q are positive. 

If v(x) is identically the positive constant r(a) on [a,b], then q 
is proportional to r and hence must satisfy the same equations. Thus 

(L+h/ [7] =O 
and 
“VA COSO + 714) Sim © = ~SiHnUSE +5/N6 ceso =0 
implies that 
- 2/4) 050 + ZG) smeo = 6. 


By hypothesis on z, and Zo (hence q) equality does not hold in all 


1 


conditions (24) and (25). We must then have, 

96) ces + 9b) sh > o 
from the second sti of (25) st P = 7/2. Hence q(b) > 0, 
a contradiction. Thus we are left with case (i), i.e., for some number 


c in (a,b), 
vie) > va) So ånd vx) >o Ir ZE. 


Now is / 
9-7 
Tih) = E ! 
Put W% (x) = r'q - rq'. Then from (24) and (26) 
(27) ~ A(LIANLG] O PNG >o 


and 


(28) 9 Cth] [A] = TPI qq? + gnn = 0 
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Solving for -qhr in (28) and substituting into (27) we see that 

/ ea x 
satisfies the differential inequality pagh E » (Lith la] 20, 
On the interval [c,b] g £ M by hypothesis, thus 


d _ V 
a in y J Y > - K DE Me 


Solving for y we may write 


M=) 
Pir) = Ae 
and for x = c 
bre) > A ) 
thus MIO) 


Vit) > Meje 


In particular 


— M (x-2) 
bib) > Wee) e E 


since v(c) A be. o i 

ath) gh) - Zov 20. 
But q(b) = 0 and q'(b) É O, therefore r(b) > O. Hence we have r 
positive on (a,b] and satisfying 

(60h) fa =O 

with 

-y íq) Coso + vaj sne =o 
and since & = 7“ /2 

nb) cost + nb) Sn hd = Wb) >O, 
Hence the function 


G= te) + RR) 


45 


satisfies the requiremencs ot Theorem 11l, and we have 
21k} = Oi) < AN) 
On the other hand, if q(a) = 0 and q(b) > O, analogously we 
may show that the solution of 
(4th) [s] = 
sth) = Cos d and sCbl= - smf 
is positive on [a,b), so that the function 
wX) È BL) 4 SIA) 
satisfies Theorem ll. 
Finally if q(a) = q(b) = 0, we have r > O on (a,b] and s > 0 
on [a,b) with r(a) = s(b) = 0; hence 
Miki VCH) FSX) 


satisfies Theorem 1l. This concludes the proof of Theorem 12. 


VI. APPLICATIONS OF APPROXIMATION FUNCTIONS 


Through the preceeding sections we have seen that the results have 
been essentially theoretical and that examples have not always been 
fruitful or plentiful. Sections 1 - 4 present one-dimensional maximum 
principles and their use in some differential equations theory while 
Section 5 actually gives some practical techniques to apply in boundary 
value problems. The problems we are considering are second-order dif- 
ferential equations with variable coefficients which in general must be 
solved by series or approximation methods. 

There does not seem to be a "best" form for an approximation func- 
tion since, in general, each boundary value problem needs individual at- 
tention. However, there is a general procedure to follow to determine 
bounds and behavior patterns of the solution (under certain conditions) 
and then the approximate value of the solution at any point in the inter- 
val under consideration: 

e Find the straight line z(x) which satisfies de boundary 
conditions as equalities and apply the operator (Lth) to it. If (Lth) 
[z] is of one sign on the interval, then z will serve as either an upper 
or lower bound as determined by the sign of (Lth) [z]. (i.e., (L+h) [z] 

£ 0 implies that u(x) & z(x) and (L+h) [z] => O implies that z(x) 

= u(x)). 

2w If the constants om 0% of the boundary conditions are 
of the same sign, then the approximation function z(x) = 0 will show 


that the solution u(x) is either non-negative or non-positive. (i.e., 


TEZ Y, > O implies u œ ORN O, UE 0 


implies u = 0). 


41 


3. If the general behavior of u(x) is desired, it is only 
necessary to bound u between the straight line from (1.) (assuming this 
is possible) and some constant function which satisfies the approxima- 
tion boundary conditions (9) and (10), or (11) and (12) in Section V. 

4, If an approximate value of the solution u(x) is desired 
at some point in the interval then find a polynomial of exponential func- 
tion that satisfies the approximation boundary conditions mentioned in 
(3.) above. If a more exact approximation is desired, then the form of 
the approximation function is more complicated and the difficulty of evalu- 
ating and checking the restrictions increases. In many instances it 
would be as easy to find an approximate or series solution (if initial 
conditions are known) and evaluate it numerically on a computer. The 
value of the above technique is that the behavior of the solution of a 
complicated second order differential equation can be determined easily 
and in less time than other methods. 

The following examples will help point out the applications of 
Section 5 and Theorem 12 specifically. The choice of the approximation 
functions may seem ambiguous, but they readily become more understandable 
in the context of each problem. 


Example 1: Consider the boundary value problem 
O We 3u'-XH =O OW (O) 
with 


o -uro+rulelj=o and w=}. 


As noted in (2.) of the preceeding discussion the solution of (1) satis- 


fying (2) is non-negative. Consider the straight line Z2/¥) = d (xt!) 
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which satisfies (2). Then 
) = M 
(Lthjl2] = 2” +32'-x2 
eh = > 
= 4 [3- x+] 2 © 
on (0,1). Hence let z(x) = Z, (Xx) and we have 
Likt) € UM) 
We note that y =. 1 satisfies the requirements for an upper bound, i.e., 
O 
(Lth) [y] & O0 and y(1) = 1, -y$ 40) + y(0) =1 2 0. We have quite 


simply determined the general behavior of u in (0,1), i.e., 


f 


(1/2) (x+1) £ u(x) £ 1. 





Fig =z 


To determine more accurately the bounds for u at any point x, we :choose 


a more complicated function for z the upper bound. A natural choice 


1? 


to consider would be 


z, (x) = (1/4) (x+3) 
which satisfies 
z(1) = 1 and -z'(0) + z(0) = 3/4 > 0. 


But z(x) does not satisfy 


(3) (L th) £2,] E é /3- 1(442))] £ (O 
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on the entire incerval (0,1). But we may define 
aa a. (tee Sy 
s [8143], odxe 2h2i-3] 
dh | : | E 
m (AT 3) ) ha aax € 4 
where (1/2)( y 21-3) is a zero of 3 = x(x+3) = 0. Then 21 is continu- 


ous and satisfies the conditions for an upper bound and we have 
1 ter) < uta) e é ja 43] 
on (0,(1/2)( 4 21 - 3))and 
dod) é ulr) € 4 (443) 
on ((1/2) ( ¥ 21 - 3,1). 
| Example 2: Consider 


(4) u" - 14x”) u ="0 tor U < x <1 


with 


(5) u(0) = 1 and u(1) = 0. 


The solution u must be non-negative. Consider the straight line 


z| (x) = 1-x 
which satisfies (5). Now 
(l+h) 2] = 2,0 — (118) 2, 
= O — (nx?) (1-4) £0 
since 0 < x < 1. We now have an approximate idea of the shape of u. 
We pick for 20 (the lower bound) a parabolic function. We choose z 


2 


such that Z, (0) = ] and z, (1) = 0, and so that z, is concave upward. 


2 
In Example 1 we showed steps 1-3 of the general method and here we show 
how a better bound can be found by choosing a more complicated approxi- 


mation function. Let 


9 
z, (x) = (x-1) 


which satisfies (5). Then 
(bah) fa) = 2; — (14x?) 22 
= 2- (4x?) (k~I)" 2 0 on (61). 


Hence 


RUDE! Send), = = 


In particular x = 1/2 


L L y/d 1 

A 4 

7 e uli) € &. 
Choosing for Z, the expotential function 







o mo 2) =i— y% 
Z(kK)= e (\-x) 


which satisfies 
2 (0)=| and 2(i)=0- 


Now 


(Lt) [3] = 2 — C4 x) 2, 


H 


= - X 
(3-k)e o (4x2) (= Ade 


o* [x®-x* +2] 20 61 (41) 


if 


since e * >0and 2 > x - x). Hênce 


EM ID io 4 


and in particular at x = 1/2 


303 £ u(x) £ .500, 


a considerable improvement. 


Since 2, £0) =z, (Ù) = u(0), we may bound u' (0). Now 


4r Iž) — lo) 
Alo) = own | EA 


) 


and from our previous results 


-2 = 2/0) € U0) € 20)=—1, 


Example 3: Consider the problem 


(6) WY =-xH =O FOR O<XK< | 
with 
(7) nervo r nual =o and Wh = / 


Noting the similarity of boundary conditions (7) with boundary condi- 
tions (2) in Example 1, we choose 
| 5i 
Zltje 2 (K+) 

and note that 

(hth) [aJ= 2" - x2 

-Sytt+)£o en (01). 
Hence Zz, (x) = (1/2) (x+1) is an upper bound for the solution u(x). We 
choose for Z, the exponential function 
2 s | 
Zl) = 
which satisfies (7). Then 
FO] 

(L+h)laJ = (xe 20 2 (01). 

Thus 


e a hix e $x), 
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and in particular at x = 1/2, 
0.6065 &u(l/2) & 0.7500. 
An interesting application of the theory developed allows us to 
bound solutions to the differential equation 


(hr alu] = wte gusu'+ ARU = hk) 


on (a,b) even when the coefficients g(x) and h(x) are piece-wise continu- 


ee 
Example 4: Consider 
o u"+xu«'— hirn =0 Im ene 
with | 
(9) wlit=o ana mino! 
where 


|: 04x £ 1/2 


h(x) = 
2.5 9182 & x“ Gel 


From the boundary conditions consider the approximation function 


z(x) = x 


which satisfies (9). Then 


(44h) [27 


dj 


P+ LP! — Afr) 2 


KTE x «os “exe 


s 


Hence 0 4 u(x) & x, since u must be non-negative by (9). A better 


lower bound may be determined by considering the parabola 


2) (x) = i 


which satisfies (9) and 
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i 
fu 
R 
x 
su 
| 
7 
X 
Ny 
po 


(Lth) lz 


Thus 
a & u(x) € x. 

When working through the theoretical background for this technique, 
there is a question of generalizing the approximation functions so that 
the "best" bound for the solution u can be obtained, In considering this, 
the following discussion will point out some problems which arise. 

We will choose general expressions containing parameters for the 
approximating functions z, and z, and attempt to pick values of the 


] 2 


and Z5 satisfy the appropriate conditions. For 


A [2- € 0) 
g J2 E eee 


where A,B and œ are constants to be determined. From Section V we 


parameters to make zy 


example choose 


$? 


2, it) 


and 


i 


a (x) 


know that z (x) must satisfy 


(4+h)[2,] < fx) on (a,b) 


and 


20)> 14 5 26/264, . 
Then 


(hth) lR] = 2 +9 10)2' +h 
E es _ 1X “ea I) AAA) pent) 


Now choose Xx so large that 

x2= of G(x) thik) > ©, 
where f,g and h are assumed to be bounded on [a,b]. To make z (x) an 
upper bound, we must have 


— X [X-A 


) 
(hihi) = 2Ahi -= A(x x9 in) hits) e < A, 
| — Linn Aj 
A [2 hiy = (K aginithinye Telm). 


This is satisfied if we let 


; — ol (Xd) 

b= mn [-2hit) + pe eg) bre 
hi keh 

and then if we choose A such that 


/ | 
A> £ ma - rot. 
OexEl 
By its definition k >0. We have now satisfied the condition 


(44h) [2,1 elx) br aex<b, 


We must also satisfy 


-K aA) 


z(a)= A/2- € = A21 O) 


2, (b) = Al 2- RJ 27%. 


m xl b-a) 
Let c=2 «e » and choose A such that 
ik / 
A= makiy da i mažfef, of 
o) CEK Sf 


i.e., A is the largest of the four numbers in the brackets. 


Then 


ter, 


u(x) & z, (x). 
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Using the same x in z, (x), we choose B as the smallest of the four 


numbers 
C gexel 
then 
— A(X-A) — x (X= 4) 
BIZ-e le U(x) & ft |2~e | 
on [a,b]. 


Example 5: We solve the following problem in two ways. Consider 


u! - xu = 0 for0 <x < 1 
with 
u(0) = O and u(l) = 1, 


Let z, (x) = x which satisfies the boundary conditions, and 


(Lh) [z] = -A L aoon ON | 
Hence 
0 £ u(x) <& x. 


For 2, we choose 


i(k) = X— Byli- x) 


which satisfies the boundary conditions, and B > 0. Then 


(44h) l2, [> -x*+8|24x*(-1)]> O (0,1) 


for B > 1/2, which can be seen by considering the values of B which 


make 


(444) L227) =O 
Az! 


(lth) [2]/ >O 
K=O 


since 
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and we are dealing with a cubic equation. Hence we choose B = 1/2, and 
we have 
(1/2) (etx?) £ u) ó xy 


since z on (0,1). From the theory developed before this problem, 


Ca 
Ra 


choose the approximation functions as follows: 


z= A fz- 7] 


20)= B612- e] 


Then we must choose cx so large that 


and 


ue so on (0,1). 


hence choose CX = 1. Then 


-x 
k = min [2% + (1- -x je i= (at x = 0), 








O£X4 | 
and 
c= 2 = a 
Hence 
| 
à = max] 0 ) O, JE o 
2 se ŽE =| 
and 
E 
B = aia [0z -; %0] = 
Hence 


ć Ux) e a E. |2-e E 


which does not improve on our previous results. 
Tne method outlined previously and illustrated in the example also 


leads to very interesting aspects of bounding solutions of second-order 


boundary value problems with general boundary conditions. 


We may find the straight line 
ao ZK} = COX% +e 


which satisfies the general boundary conditions: 
— 2'lajwse + 244) 3hne2% N 


(11) 
24o) sP + Hb Snt IN 


if the parameters © and + satisfy certain restrictions. 


We may rewrite (11) as 
0, (-6056) + dard) sine y 


FL 


il 


Ca COSA + (0b te) Sin 


or equivalently 


à jasme - 256] + 4 sin & 


i 


% 


a f bsn g +005 Pb | +A Im br ha, 
which has a non-trival solution for c, and c, if the determinant of co- 


1 2 
efficients 


De = 


(A Sine —Cosg) sme 


(b singh #0036) Siva 
is non-zero (assuming ME + ee Æ 0). 
D= (a~ b) sine sin A - [sin @ Cos fe + cose sing / 
= (a-b) IM o sing “sin (oe+7) o 


Hence if & and F are not solutions of 


a2) Sin(6+#) = (a-b) Sina sm 
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then we may always find the line (10). We note that if © + A = ny” 
and O =nD , f = ni“ for (n = O, 1,...)ž then (12) is 
satisfied. 

In Example 1 we used a polygonal line as an upper bound. This 
procedure may not be fruitful in the non-homogeneous case. The reason 
for this will be made clear later. 

Example 6: Consider the non-homogeneous boundary value problem 

nt = CRE 
with 


HO)= UC) = o. 


To determine possible concavity of the solution, consider 
se E 


then 


(lth) [ele 2"- B= -/ < 270-4" en (6,1) 5 


— 
= 
— 


hence z 1 is an upper bound and the solution u(x) is concave upward 


since 
Pe aa on (0). 


To determine a lower bound, consider 


MN) = 


then 

(44) [2] = 22> 3 >24X-x on LW); 
hence z = -3 is a lower bound since the appropriate endpoint in- 
equalities hold. 

By substitution we may verify that 


Upjl= #*- x% 


is the solution. We might expect then that the polygonal line 


39 


—-Xi, ogre 
21%) = 4 —h) Geox oy 


would satisfy the theoretical conditions for a lower bound since 


z, (x) & u(x) on (0,1). 





Fig. 4 


Note that this is the opposite of the procedure in our theory where we 
determine bounding functions by applying the restrictions to them 
rather than determining whether actual bounds of the solution satisfy 


the restrictions. Recall that for 2, to be a lower bound, then 


(i) (L+) [22] > f(x) = 2+x - = on (0,1) 
(ii) z,(0) < 0 and 2, (1) < 0. 


Applying the operator (Ith) to Zs 


(Lh) (2) = 2, — 2. = K’ OLK è 4 
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but 


2 (0) = Ron Lie 


From our theory this would seem to imply that z, is an upper bound for 


Z 
u; but by an oversight we have not considered the fact that each segment 
of the polygonal line must satisfy the restrictions individually. 
Hence for z, (x) = -x [0,1/4] we see that 
(Lth) l] = x< fixy 
with 
2,/0 =0 Amd 22(%) = -h eo. 
Thus our theory does not apply. We find similar results for the re- 
maining segments of Zo. 
This points out an important result that for the non-homogeneous 
case the theory developed may not apply to the polygonal line approxi- 
mation functions, 
The reason for this is clear upon consideration of Example 1, where 
we choose the straight line 
24) = 2 (143) 
as an upper bound, and find that 
(lihlla la 2"132'-x2 = É /3-x(x+43)/ 
changes sign on (0,1), i.e., (Lth) [z] has a zero in (0,1) and by 
suitably changing z(x) we can satisfy the restrictions 
(44h) [2] €0 on (G1). 
But for the non-homogeneous case, we must satisfy 
(Cutb)fe) < Hx) 
and we may not be able ro suitably change our approximation function 
to satisfy the conditions. 
This points out the importance of parametric approximation func- 
tions and being able to choose the parameters in such a manner as to 


satisfy the restrictions. 


Consider then 
me O 
am=A[2- € | 
~o X 
2X = B [2 = 
and choose A, B and& according to the theory developed preceeding 
Example 5 so that Z1 and z, are upper and lower bounds, respectively. 
Then o must satisfy =“ 2 4 1 > 0; hence, choose X = 2, Then 
choose 
2X 2 
k= mn |z+ 3e = 24 3e at k=1. 


04% E/ 
We must now determine 
max -twf = Mex gera fo -7 af Xx=0 m~ a 


ospel vt xe 


and 


mer ff = mom eye be e at x= % , 


OENE] 


Since €, = 0, we do not need c = 2 =e since Bote = 0. Thus 


choose A and B so that 


= Zz 
A = mp | 6 O, = 2 | = 
2732 














and 
B= min /20, =» O ja - i 
| Gere * ? grizo 
Thus 
= |e] Z E, 
G+12€7* 


We may graph the above by approximating the function Z, (x), and we 


have o | 
HX) 
ei 
2,(1) 
-2 
Pip go 


We now have the option of attempting to change the parameters A, 

B and “ to obtain better bounding functions 2, and 25 or keeping the 
existing functions and tolerating the possible large error. In either 
case we are sacrificing time or exactness and are defeating our pur- 
pose. 

The technique used in the preceeding examples has indicated the 
possibility of approximating solutions u(x) by multiplying the straight 
line s(x) found above by "warping" functions & (x) and S (x) so that 

PAS é nae J S 
This technique should not oniy help determine the general behavior of 


the solution but yield usable numerical approximations. 
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